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Abstract. The elastic and capillary interactions between a pair of colloidal particles trapped on top of 
a nematic film are studied theoretically for large separations d. The elastic interaction is repulsive and 
of quadrupolar type, varying as d~^. For macroscopically thick films, the capillary interaction is likewise 
repulsive and proportional to d~'' as a consequence of mechanical isolation of the system comprised of the 
colloids and the interface. A finite film thickness introduces a nonvanishing force on the system (exerted 
by the substrate supporting the film) leading to logarithmically varying capillary attractions. However, 
their strength turns out to be too small to be of importance for the recently observed pattern formation 
of colloidal droplets on nematic films. 



PACS. 82.70.Dd Colloids 



i.OS.Cd Surface tension and related phenomena - 61.30.-v Liquid Crystals 



1 Introduction 



The interactions of colloidal particles trapped at fluid in- 
terfaces have been found to differ significantly from the 
corresponding interactions in bulk solvents. This has been 
studied mostly for electrically charged particles trapped 
at interfaces with water. On one hand, the presence of 
the interface gives rise to direct dipolar electrostatic re- 
pulsions between the colloids (see Refs. [I][2ji3j for some 
experimental evidence), on the other hand deformations 
of the interface may induce longer-ranged capillary at- 
tractions (briefly reviewed in Refs. [UllllSj) which is pos- 
sibly the source of pattern formation observed in various 
experiments [TllSligifTOlITT] . (See, however, Ref. [H] for an 
alternative explanation due to interface impurities.) 

Recently [T31IT4]. the experimental observation of or- 
dered structures of glycerol droplets bound to a nematic- 
air interface has been reported and attributed to an ef- 
fective pair potential between the colloids which contains 
a repulsive, elastic part due to director deformations in 
the supporting nematic film and an attractive, capillary 
part which is long-ranged and mediated by logarithmi- 
cally varying deformations of the nematic-isotropic inter- 
face caused by the droplets. The schematic setup of this 
experiment is depicted in Fig. [TJ According to Ref. [T^ . 
the colloidal particles experience an upward force caused 
by elastic forces due to director deformations in the sup- 
porting nematic film. This upward force on the colloids 
is supposed to give rise to the aforementioned logarithmi- 
cally varying interface deformation. Applying a superposi- 
tion approximation for the deformation field, one can show 
that the ensuing effective capillary interaction potential 



between two colloids is likewise varying logarithmically. 
This is similar to the flotation interaction of mm-sized 
particles at fluid interfaces for which the force on the col- 
loids is caused by gravity (see, e.g., Ref. [15]) and also 
parallels the tentative explanation given for the exper- 
imentally observed attractions between sub-^m charged 
colloids at a water-oil interface [16] (for the controversy 
around this explanation see Refs. [T7[[TSl[Tgi[^[^Tll^^ll^ 

m)- 

However, it is now well established pTllTHlfTM^Tll^ 
[^[^ that interface deformations and effective colloidal 
interactions varying logarithmically only arise in experi- 
mental systems which are not isolated mechanically. For 
mechanically isolated systems it can be shown pSl[?7l[^ 
[^ that both the interface deformation around a single 
colloid and the effective interaction between two of them 
are shorter-ranged and the latter cannot be calculated re- 
liably within the superposition approximation. 

In the following we will extend the arguments pre- 
sented in Refs. ^1^1^ 7.28,29 to systems with colloids 
at nematic interfaces. We will show that mechanical iso- 
lation of the system "nematic film - colloid - air" can 
be violated through a subtle interplay between the finite 
thickness of the film and the anchoring conditions at the 
colloids and at the nematic interfaces with the substrate 
and with the air, respectively. However, for experimental 
conditions as the ones described in Ref. [13] a quantita- 
tive estimate of the strength of the ensuing logarithmic at- 
traction between the colloids yields that these attractions 
are unobservably small. Therefore it seems likely that this 
kind of asymptotic capillary forces cannot be invoked as 



M. Oettel, A. Dominguez, M. Tasinkevych, S. Dietrich: Effective interactions of colloids on nematic films 



m 



air 




nematic film 



A A A A A A^ 

\ \ \ \ \ V ^^^^v- 



glycerol substrate 



V A A A A A A A A 

^ \ \ \ \ \ \ \ \ 



Fig. 1. Scliematic setup of the experiment reported in 
Ref. [13]. Colloidal glycerol drops {R — 1 ... 7 ^m) are trapped 
at the surface of a thick nematic film {h ?a 60 fim). The di- 
rector field in the nematic film is sketched by the black lines 
and their dotted interpolations. Nematic anchoring at the glyc- 
erol substrate at the bottom of the film as well as on the sur- 
face of the glycerol drops is parallel, whereas it is perpendicu- 
lar (homeotropic) at the deformed nematic-air interface. Each 
drop is necessarily accompanied by a topological defect ((g)). 



a relevant mechanism to account for the observations re- 
ported in Ref. [13] . 

The manuscript is organized as follows: In Sec. [2] the 
coarse-grained model for the nematic phase is introduced 
which will serve as the basis for all subsequent calcula- 
tions. In Subsec. 12.11 we study the case of an infinitely 
thick nematic film. First, we compute the asymptotic form 
of the director field and the ensuing elastic force between 
two particles. Then we calculate the effective force aris- 
ing from the deformation of the fluid-nematic interface 
caused by the clastic stresses. In Subsec. [2?2] we consider a 
nematic film of finite thickness, which models more closely 
the experimental setup described in Ref. |13| : for such a 
system we extend the above calculations to the two oppo- 
site cases of perpendicular and parallel anchoring of the 
director field at the substrate surface. In Sec. [3] we discuss 
our results. 



2 Coarse-grained model 



In view of the mesoscopic length scales involved we de- 
scribe the bulk part of the nematic free energy associated 
with the director deformations in terms of the Frank free 
energy expression within the one-coupling approximation 



•^nc = / d\ f\r) 

/ d^r [(V • nf + (V X n)2] 



K 

— / cPr Wui ■ Vti, 
K 



(1) 



/ dV V • [n(V • n) - (n • V)n] 



V^nc denotes the volume occupied by the nematic film, n 
is the director field (n'^ = 1), and the constant K is of the 
order of 10"^^ N [30]. The total divergence term in the 
last line of Eq. ([1]) exhibits the so-called "if24~structurc" 
and is unimportant for the bulk equations describing the 
equilibrium configuration. The surface free energies asso- 
ciated with the interfaces with air and substrate, respec- 
tively, (see Fig. |T]) are described in terms of the Poulini 
expression |31| : 



K.-^ dA(n.e^)2 + ^ 



dA (n • ba) 

^sub-no U -^coll 



(2) 



Here, e^i denotes the local surface normal unit vector 
pointing outwards from the film or the colloid. Normal 
alignment is favored for Wi < and parallel alignment 
for Wi > 0. Typically one has \Wi\ - 10"^ N/m [31] so 
that the length scale K/\Wi\ ~ 1 ^m is smaller than the 
range of droplet radii investigated in Ref. [13j . Thus, in the 
"strong anchoring" limit which we shall consider, the effect 
of the boundary terms is so strong that as a first approxi- 
mation it amounts to fixing the angle between the director 
and the surface normal. We shall adopt Wi < (normal 
alignment at the nematic-air interface) and W2 > (par- 
allel alignment at the nematic-glycerol interfaces). Some 
consequences of deviations from the strong anchoring limit 
will be discussed in App. [XI 

These surface contributions to the free energy ("wet- 
ting energies" ) are small corrections to the surface tensions 
which are mainly due to dispersion interactions. We de- 
note these non-nematic contributions to the surface ten- 
sions as 71 (colloid-air surface tension), 72 (substrate- 
nematic surface tension) and 7' (nematic-air surface ten- 
sion). Typically, these surface tensions are of the order 
of 10~^ N/m. Therefore, they are much larger than the 
constants Wi which determine the nematic contributions 
to the surface tension. We see that due to the above an- 
choring conditions the surface tension of the substrate- 
ncmatic interface carries no nematic contributions due to 
n • e^ = (Eq. ^) whereas in the strong anchoring limit 
the full nematic-air surface tension is 7 = 7' -f Wi /2 w 7' 
due ton- CA = 1 (Eq. ^)^ 



^ A genuine contribution to the surface tension of a nematic 
interface arises if one takes into account the variation of the 
nematic tensorial order parameter through the interface, de- 
scribed by, e.g., the Landau-de Gennes free energy functional 
(generalizing Eq. ([T}). The magnitude of these contributions 
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The canonical stress tensor TTy associated with the free 



energy expression in Eq. ([T]) is given by 
9/b 



9nfc, 



nk,2 - 5ij f 



(3) 



where Uk^i = duk/dxi (summation over k). The total 
stress tensor 77^ is obtained by adding the contribution 
of the isotropic pressure p: 



n, 



S^J P 



(4) 



2.1 Macroscopically thick nematic film 



First we consider the limiting case h ^ oo (i.e., very thick 
nematic films (see Fig.[T]))oDue to the small values of the 
elastic coupling constant, K ^ ^ R, and of the anchoring 
energy, \Wi\ <C 7, the equilibrium configuration of a single 
colloid at the nematic-air interface deviates only slightly 
from the reference configuration depicted in Fig. [21 In this 
latter configuration, the interface is flat and the colloid is 
positioned such that the contact angle fulfills Young's law 
cos 6* = (71 -72)/7- 

The total force on the whole system reads (the super- 
script '^^^' denotes evaluation on the positive (negative) 
side of the oriented surface, i.e., on the side the arrows in 
Fig. [2] [which indicate the surface normals] point to (do 
not point to)): 

fdA-n^f dv{w-n)+f dA-{u+-n-) 

Js Jvi u V2 Js„,^„ u Si 

= I dA (tt^^ + Pair -p)e^ + 



dA-[-K + (pair -P)l] 



Si 



dA TT^^ e^ + / dA- 



(5) 



In obtaining this equation we have applied Gauss' theo- 
rem. Furthermore we have used the relation V • 11 = 
in volumina Vi and V2 which is valid because the ref- 
erence configuration is taken to be in force equilibrium. 
This also implies that the isotropic pressures above the 
interface (pair) and below it [p) are equal and that the di- 
rector configuration is given by the corresponding Euler- 
Lagrange equilibrium equations following from the func- 
tional in Eq. ([T|). Since at the interface Si the colloidal 



can be estimated by the surface tension of the interface be- 
tween the nematic and the isotropic phase of a liquid crystal. 
Typically, such a surface tension is also 0{1Q~'^ N/m) ~ \W'i\ 
and therefore small compared to the dispersion force contri- 
bution. Furthermore, a distorted director structure in the bulk 
may also give rise to surface energy contributions on the bound- 
aries. For an example see Ref. [30], p. 131 and p. 174. Also in 
this case it can be argued that the corresponding contributions 
to the surface tension do not exceed \Wi\. 

^ This was implicitly assumed also by the authors of Ref. [13] 
in discussing Fig. 2 therein. 




nematic phase 



Fig. 2. In the reference configuration the whole system is 
divided into volumes Vi and V2. Volume V2 (enclosed by the 
upper dashed curve) includes the air and the glycerol drop and 
volume Vi (enclosed by the lower dashed curve) includes the 
nematic. The arrows indicate the direction in which the sur- 
faces (including the infinitesimally displaced ones) are oriented: 
S encloses the whole system, 5mon is the interface between the 
nematic phase and air (acting as a meniscus), and Si is the in- 
terface between the colloidal drop and the nematic phase. The 
director field and the topological defect ((g)) are indicated as in 
Fig. [1] The radius of the three-phase contact line is denoted 
by po = -R sin 6 where 6 is the contact angle of the air-nematic 
interface with the colloid of radius R. p denotes the lateral 
distance from the vertical symmetry a^ds of the colloidal drop. 



drop is rigidly attached to the liquid crystal, we can iden- 
tify the vertical force F on the colloid and the total force 
Ft^ on the air-nematic interface by 



and 



Fe, 



F^^- 



dA- 



Si 



dA TT, 



(6) 



(7) 



respectively. Mechanical isolation of the system means 
that the total force ^„ dA • 11 acting on it is zero which 
leads to 



F^F, 



(8) 



For a given force on the colloid and a given stress on 
the interface, the interface deformation relative to the ref- 
erence configuration can be determined perturbatively To 
that end, we summarize briefly those results of Ref. [21] 
which are pertinent also for the present system. With the 
introduction of the two small, dimensionless parameters 



EF 



-F/{2ttpoj) and 



-i^„/(27rpo7) (9) 
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one can expand (up to second order in ep and Ett) the 
free energy difference JF associated with the interface de- 
formation u{p > po) around a single colloid (see Fig. [2]) 
and with a vertical shift Ah which is the difference of 
the colloid center position relative to that in the reference 
configuration: 



T. 



27r7 / dp p 

'pa 



du\ u^ 

Tp) ^2}? 



1 



7r7[u(po) - Ahf - FAh 



(10) 



Here, A = {'^/{gpm)Y^'^ is the capillary length associated 
with the interface where g is the gravitational constant 
and pm is the mass density of the nematic phase. This 
expression for the free energy contains all surface free en- 
ergy changes relative to the reference configuration involv- 
ing the interfaces between air, nematic, or colloid. It also 
contains the contributions due to volume forces acting on 
the nematic (associated with A) and the energy change 
of the colloid upon vertical shifts (for further details see 
Ref. [21] ) . Note that to leading (quadratic) order in e^ , ep 
the free energy change of the nematic due to the shifted 
interface and due to a change in the director configuration 
with respect to the reference configuration is captured by 
the term oc / tt^z u. (The analogous textbook argument for 
electrostatics [35] can be easily generalized to the nematic 
case described by the free energy expression in Eq. H]).) 

Minimizing J- with respect to u{p) and Ah and fo- 
cussing on the regime p ^ A yields 

w(p) - Po(e7r -eF)ln V- / do- cr7r22(cr)ln- ,(11) 

P 1 Jp P 

with C ~ 1.12. We see that in the case of an isolated 
system (e^ = ep) the logarithmic part of u{p) vanishes. 
The second term on the rhs of Eq. pT|) leads to u{p -^ 
oo) ex /9^"+^ if TTzz oc /9^" and thus describes a shorter- 
ranged power-law decay of the interface deformation. 

The absence of logarithmic deformations for an iso- 
lated system has been derived here under certain simplify- 
ing conditions (small interfacial deformation everywhere, 
rotational symmetry). In App.lBJwe demonstrate that this 
conclusion holds in general. 



2.1.1 Asymptotic director configuration and elastic force 
between colloids 

In Refs. [551IM] it has been shown that a colloidal drop im- 
mersed in the bulk of a liquid crystal is accompanied by a 
single counterdefect such that the total topological charge 
is zero (here, the volume occupied by the colloid contains a 
topological charge which may be represented by a virtual 
defect inside the colloid) and the asymptotic behavior of 
the director field is of dipolar character. Based on similar 
considerations we shall show that for a colloidal drop lo- 
cated at the air~nematic interface the boundary conditions 
for that interface impose a quadrupole-like asymptotic be- 
havior of the director field. Macroscopically far from the 



colloid the director is oriented parallel to the z axis. Ac- 
cordingly, at large but finite distances r the director is 
given by n(r) ~ {ni,n2, 1 — 0(71.^,^2)) and the bulk free 
energy corresponding to Eq. ([1]) is given by 



TTD ^^ 2— 



d^r 



V„c 



J2 iVn^f + Oinj) 



(12) 



Here we have discarded the total divergence term in the 
free energy expression ([IJ. It is unimportant for the bulk 
equations and it adds a mere constant to the free energy 
because the director is anchored normally at the bound- 
ary (the nematic-air interface). Thus for each component 
i = 1,2 the equilibrium director field fulfills the Laplace 
equation 



Ar 







(13) 



Analogously to electrostatics, the asymptotic solution for 
Ui can be expanded in terms of multipolcs. To this end we 
consider the reference configuration in Fig.[2| We choose as 
the origin of the coordinate system the center of the circle 
formed by the planar three-phase contact line. The solu- 
tion for the director field has to fulfill the following require- 
ments: (i) rotational covariance around the z-axiqj and 
(m) ni{x,y,z = 0) = 0. Analyzing the multipole ansatz 
(with r = (ri,r2,r3)) 



1 + J2 P^c ^ 



tQ../-^ + ...iU) 



it follows that rotational covariance requires qt ~ 0, Pia = 

-I OtQ + ^mag £iq3 and Qiap = Q p Oia 4- Qmag^p ^ia3 i^ijk 

is the Levi-Civita tensor). Pmag and (5^ag/3 ^^"^ dipole 
and quadrupole moments, respectively, for a director field 
of "magnetic" type, i.e., for which divri^ = holds. The 
boundary condition (ii) at the interface with the air fur- 
ther imposes P = P^ag = and Q^ = Q Sps, Qmag,/3 = 
Qmag ^/33- It appears to be difficult geometrically to match 
the asymptotic solution of "magnetic" type with a solu- 
tion near the colloid which obeys parallel anchoring at 
the colloid surface. Therefore we discard the magnetic 
quadrupole, i.e., the leading asymptotic term is given by 
the remaining, "electric" quadrupole term: 



zr\ 



(z = rs) 



(15) 



This is at variance with Ref. [13] (see Eq. (3) therein 
and the considerations in the paragraph above that equa- 
tion which assume a dipole field) but it is consistent with 
the analysis in Ref. [34] . Dimensional analysis yields Q = 
0{R^) [31]. Note that we have derived the asymptotic be- 
havior of the director field at the interface using strong 
anchoring at the interface (rii = 712 = 0). In App. lAl we 
discuss corrections to strong anchoring which are, how- 
ever, of subleading character and leave the leading behav- 
ior (Eq. P^ ) unchanged. 



•^ If D specifies the transformation matrix for such a rotation 
then this requirement is given by n(D • r) = D ■ n(r). 
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The asymptotic elastic interaction between two col- 
loids in the bulk at distance d accompanied by a quadrupo- 
lar director deformation has been analyzed in Ref. [35| 
(for weak anchoring) and in Refs. [36JI34] (using a coarse- 
graining method, applicable also for strong anchoring) 
yielding identical results. For the present configuration 
(distance vector perpendicular to the asymptotic director) 
the elastic potential is repulsive and varies as 






(16) 



We have used that the diniensionless force parameter ep 
is proportional to Q^ which actually follows from Eqs. ©, 
(O, (O, and pT|) below. Note that we have simply extrap- 
olated the bulk results for two interacting colloids which 
cause asymptotically quadrupolar deformations of the di- 
rector. This appears to be reasonable because the asymp- 
totic director field in the nematic phase for the interface 
problem is precisely that of the bulk solution in the lower 
half plane and because the bulk solution is antisymmetric 
with respect to z -^ —z, thus respecting the boundary 
condition ni{x, y,z = 0) = Q. However, the precise numer- 
ical value of the quadrupole moment Q might be rather 
different for the case of the colloid trapped at the interface 
as compared to the bulk case. 



2.1.2 Asymptotic behavior of the stress on the interface and 
meniscus-induced effective potential between colloids 

The asymptotic behavior of the stress tensor component 
-Kzz at the interface follows from inserting Eq. (fT51) into 



Eq. ©: 



^^ I interface 



fE(. 

2 ^ p8 



2 _ 2 ^ 



z=0 



(17) 



lines dSi. The meniscus-induced effective potential is the 
difference between the equilibrium free energy of the two 
colloids at distance d and their free energy at macroscopic 
distance: 



V^non(rf) = ^oq(rf) ~ ^cqid ^ Oo) 



(19) 



As before, minimization with respect to u{p) and Ahi 
renders the equilibrium free energy. 

The behavior of Vmcn{d) has been analyzed in detail 
in Refs. [21,29 . Here we summarize these results as far as 
they are relevant for the present problem. The interfacial 
stress TTzz may be decomposed generally as 



TTzzip) = n,,{\p - p^\) + n.^ilp - P2\ 

= T^zzA + TTzzS + 2 -Kzz.ui ■ 



■2n, 



.(P) 



(20) 



Here, TTzz,i denotes the stress around colloid i which per- 
tains to the problem of a single colloid. To quadratic order 
the asymptotic director field around two colloids is given 
by the superposition of the components n^ of the single- 
colloid solutions and thus to this order we recover the 
decomposition in Eq. (|20p with the mixed component of 
stress field Trzz,m given by 



_ :^ ^2 (P-Pi)- (P-P2) 



P1\'\P- 



P2I 



(21) 



It turns out that for a system under an external force 
(stt y^ sf) the mixed term TTzz.m does not contribute to the 
leading term in Vmon- Thus in the case of a non- vanishing 
external force this leading contribution to Vmon is obtained 
by a superposition ansatz which consists in approximating 
the interfacial deformation and the total stress field by 
the sum of the respective single-colloid quantities only 

{tTzz ~ TTzz.l + TTzz.li " » "l + "2) [Slllin] : 



Kicn(po < d < A) ~ -27r 7 po (e^ 

(e^ ^ £f) ■ 



epfln^ (22) 



ip' 



rl). Consequently the interface deformation Yot an isolated system (e^ = ep) and for the stress given 



around a single colloid for a mechanically isolated system j^ gq, ^ the superposition approximation tt. 



obeys u{p — > oo) a p~^ (see Eq. pT|) ) 

For the problem of two identical colloids located at p^ 
and P2 (vectors arc defined in the interface plane z = 0) 
a distance d= \Pi— p2\ ^ R apart the expression for the 
free energy is a straightforward generalization of the one 
for the single-colloid free energy given in Eq. (fTO|) : 



T^zz,l + 

T^zz,2 yields Vmon oc £p/d^ as the dominant term. How- 
ever, this is not the leading term, which rather stems from 
TTzz.m- This term has two peaks around the colloid centers 
and therefore close to the colloids it can be approximated 

by 



'S,„e„ 



\yu\ 



u 
2A2 



(18) 



2d^ 



E(-i) 



i=l 



edjp- P^ 

\p-p.\' 



(23) 



E 

1=1,2 



7 



^PO JdS., 



de \Ah, 



F,Ah, 



Here, Fi denotes the force on colloid i and Ahi is the rel- 
ative position of its center. The integration domain S'mcn 
is the whole interface plane except for the two circular 
disks bordered by the (reference configuration) contact 



where e^ = (p2 — Pi)/d- As discussed in Ref. [29], the 
qualitative behavior of Kncn('^) is captured by the integral 
over TTzz,m- 

Kicn(rf) OC -IpIef I d^pTTzz,yn{p) OC 7Po 4 (^) (24) 
(Ett ^Ef) ■ 
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We note that this form of the mixed stress and thus the '^-^ 
leading behavior of Vmcn{d) is formally analogous to the 
contributions of the electric field components parallel to 
the interface in the case of charged colloids [29] . The meniscus 
induced potential T^ncn is repulsive and falls off asymptot- 
icalljQ ex d^^, as does the likewise repulsive elastic poten- 



tial V^i. We note that Vd oc lepl and V^i, 



oc e 



_F' 



SO that 



the meniscus-induced potential is small compared to the 
elastic one for the parameters of the experiment reported 
in Ref. [13]. 

In order to explain the experimentally observed attrac- 
tions, in Ref. [13| a perturbative picture similar to the one 
presented above was suggested, but the effect of the inter- 
face stress TTzz was neglected completely. (A similar error 
has been made in Ref. [16] which the authors of Ref. [13] 
refer to.) In a heuristic way, only an "upward" force on the 
colloid (perpendicular to the interface) associated with the 
anchoring "wetting" energy at the nematic-particle inter- 
face has been invoked in Ref. [13], neglecting the force 
on the interface described by tt^z- In this way, the unbal- 
ance of the force gives rise to a logarithmic term in the 
meniscus deformation and in the meniscus-induced po- 
tential. However, mechanical isolation (i.e., force balance) 
renders the meniscus-induced potential actually repulsive 
and shorter-ranged (see Eq. (jM]) ). apparently in contrast 
to the experimental results. A logarithmically varying po- 
tential can only arise if mechanical isolation is violated 
(see App. |B|. Below we shall investigate whether a net 
force on the system "colloid and interface" may appear if 
the thickness of the nematic phase is finite, as it is the 
case in the actual experiment. 



2.2 Finite thickness of the nematic film 



In our discussion of a finite film thickness of the nematic 
phase we shall consider two cases: 

1. The anchoring of the nematic director at the surface of 
the bottom substrate is perpendicular as it is the case 
at the upper interface with the air. This case bears a 
strong formal resemblance to charged colloids on water 
surfaces which have been discussed in Refs. [^T[|29j . 

2. The anchoring at the bottom substrate surface is par- 
allel (as in the experiment reported in Ref. [13]). At 
large lateral distances from the colloids, this leads to a 
director field which gradually rotates from parallel ori- 
entation at the bottom substrate to the perpendicular 
orientation at the upper interface. 

For both cases the total force on the system - compris- 
ing air, nematic film, colloid, and the substrate - must 
be zero. This leads to (see Fig. ^a.) and compare with 



'^ Superficially one would expect a leading decay Knen oc d^"^ 
as displayed in Eq. (|23|l . However, due to the geometric fac- 
tor in the numerator of Eq. 
vanishes upon integration. 
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Fig. 3. The reference configuration for a nematic film (a) 
differs from the reference configuration shown in Fig. [5] by 
the addition of the substrate volume V2- The interface be- 
tween the substrate and the nematic film is denoted by Ssu_h- 
Panel (b) shows the string of image quadrupoles Qt ~ Q which 
are needed to fulfill the boundary conditions for perpendicu- 
lar anchoring of the director field (originating from the colloid 
quadrupole Q) at both interfaces confining the nematic film. 
The distance between any two nearest neighbor quadrupoles is 
2h. 



Eq. ®) 
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23)1 . this apparent leading order 



For reasons which will become clear in the discussion of 
the second anchoring case, we consider the isotropic pres- 
sures in the substrate and air, psub and pair, respectively, 
not necessarily to be equal to the pressure p in the ne- 
matic film. The second equation in Eq. ([^5)1 follows from 
the equilibrium condition V • 11 = which holds in all 
volumina. 



2.2.1 Perpendicular anchoring at both interfaces 

As discussed above, the presence of the colloid asymptot- 
ically generates a quadrupolar director field which fulfills 
the boundary condition at the nematic-air interface. In 
order to fulfill the boundary condition at the substrate- 
nematic interface, an image quadrupole of the same strength 
Q is needed which, however, leads to a violation of the 
nematic-air interface boundary condition and requires a 
second image quadrupole etc. Continuation of this pro- 
cess leads to a string of image quadrupoles as depicted 
in Fig.|3Jb). This string of quadrupoles generates a stress 
field TT which vanishes for large lateral distances. There- 
fore the isotropic pressures must be equal in all volumina: 
Psub ~ Pair = P- From Eq. (j25p one finds that the differ- 
ence between the force F on the colloid and the integrated 
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stress F^ at the neniatic-air interfaee, 



F^ = AF 



dA -K. 



(26) 



is given by the integrated stress over the substrate sur- 
face, i.e., the total force on all quadrupoles above the sub- 
strate surface exerted by the image quadrupoles in the 
substrate. Expressed in terms of the force Fq^q between 
two quadrupoles at distance 2h, we find 

n—1 n—2 

OO ^ 

= E ;i5 = C(5) « 1.04 . (27) 

ra=l 

In this equation, the first sum is the total force (divided 
by Fq-q) on the first quadrupole above the substrate 
exerted by all quadrupoles in the substrate, the second 
sum is the total force (divided by Fq-q) on the second 
quadrupole above the substrate etc. For the force between 
two quadrupoles at distance 2h we find, using Eqs. (fT5|) . 
([3]), and 



Fn 



5 K^' 



(28) 



Thus mechanical isolation for the system "colloid and 
nematic-air interface" is violated and the total force AF 
on this system is (up to the factor 1.04) given by the 
quadrupole force Fq-q. Nevertheless the magnitude of 
the corresponding induced logarithmic capillary potential 
(see Eq. (f^ with e^ — ep = AF/(2tt^po)) is small, be- 
cause Q ^ R^ due to dimensional arguments [M]. For pa- 
rameters similar to the ones appearing in Ref. [13] {R/h ~ 
10-\ K/{-fR) ~ 10-4) we find T/mon ^ 10"" fcsT \n{R/d), 
which is unimportant for the actual intercoUoidal interac- 
tion. 

The elastic potential V^i(d) between two colloids is the 
interaction between the second quadrupole and the first 
quadrupole together with its string of image quadrupoles. 
Using the solution given in Refs. [5M55] we have checked 
that Vci{d) remains repulsive. For d < h the overall mag- 
nitude of Vei{d) is somewhat weakened, whereas for d^ h 
a crossover to Vei{d) ex cxp(— d//i) is observeqj- 



2.2.2 Parallel anchoring at the bottom substrate 

We assume that the substrate induces a preferred in-plane 
axis for the director orientation which we take to be the x~ 
axis. With no colloid present at the nematic-air interface, 
the equilibrium director field is given by 



no 



^sin(-goz)^ 


.cos{-qoz). 



qo = TT/{2h) , (29) 



^ Tliis result can be obtained more easily by solving the field 
equations Arii = in cylindrical coordinates rather than by 
using the image quadrupoles. 



with the consequence that both at the nematic-air and at 
the nematic-substratc interface a constant stress is acting: 



TTO,, 



K 2 



(30) 



For the unperturbed interface to be in equilibrium, the 
air and substrate pressures differ from the pressure in the 
liquid crystal: p - p^ub [air] = t^o,zz- 

We now introduce a single colloid at the nematic-air 
interface in the reference configuration. Due to this pres- 
sure difference, the force on the colloid and the integrated 
stress over the nematic-air interface are given by (see 
Eq. (ESD) 



Fe^ ^ dA-{n- ttq^^, 1) , 

JSi 

dA (tt^^ - 7ro,2z) 



F^ = 



(31) 
(32) 



and the total excess force on the system "colloid and 
nematic-air interface" is determined by 



F- F.^ = AF = - 



dA {-Kz2 - 7ro,22 



(33) 



In order to calculate the director field n and the stress ten- 
sor TT in the presence of the colloid, we introduce the aux- 
iliary director deformation fields v(x,y,z) and w{x,y,z) 
which parametrize the deviations from the unperturbed 
director field no and which are small at large distances 
from the colloid: 



n 



sin{—qoz + v) cos w\ 

sin w « no + 

I cos(— (7o2 -I- v) cos wj 

^cos{qoz) V + ^ sin{qoz){v'^ + w^)^ 

w 
I sin((7oz) V ~ ^ cos{qQz){v'^ + uP') . 



(34) 



0{{v,wf) 



The first equality in Eq. (|34p is a general parametriza- 
tion of the director field n in terms of the auxiliary fields 
v^w which fulfills n^ = 1. These auxiliary fields are taken 
to be zero at the substrate and the nematic-air interface, 
i.e., the boundary conditions are v{x, y, 0) = v{x, y, —h) = 
w{x,y,0) = w(x,y,—h) = 0. The nematic free energy 
of the film up to order 0(v'^,'uP') is obtained by insert- 



ing Eq. (|34p into Eq. ^ for the Frank free energy after 
dropping the total divergence of the isr24-type. Using the 
boundary conditions for v and w we obtain [37] (with the 
notation introduced in Eq. ([3|): 



T- 



■fill) 



^n + ^ 



Vfii 



d^r (wj^j + w^ 



2 2\ 



(35) 



Here, Tq is the free energy of the film without colloid. 
At first glance it is not evident that this free energy is 
positive definite. However, the boundary conditions on v 
and w ensure positivity [37j . Upon minimization we find 
Av ~ 0, i.e., the deformation field v of the director can 
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again be expanded in terms of electrostatic multipoles 
(Eq. (fT4|) ). On the other hand, the solution for w must 
fulfill the Hclmholtz equation (Z\ + (?o)w = and can be 
expanded in terms of multipoles as follows: 



qor 3> 1). On the other hand, for large film thicknesses 
there is an intermediate asymptotic regime R <^ r <^ h: 



w{qoR < qor < 1) 



w{r,B,(t>) = 



1 



E E y.^ 



(36) 



V j=Q m=-i 



v = P,- 



zx 



(37) 



The dipole contribution should vanish for /i ^ oo. If one 
assumes a power-law dependence on ft,, dimensional anal- 
ysis for Py leads to 



P.,, = 0{R^ [R/hT) (with K > 0) 



(38) 



The precise functional form of P^, turns out to be unim- 
portant for the subsequent calculations. We note that an 
asymptotic solution with a nonvanishing a;-component of 
the dipole moment cannot appear because it would not ful- 
fill the boundary conditions and the reflection symmetry 
w{x,y,z) = —w{x,—y,z) excludes any dipolar contribu- 
tion in the solution for wiQ Therefore the leading asymp- 
totic behavior of w takes the form 



{Qi J5/2{qor) + Ql Y^/2{qor)) 



+w- 



linage 



(39) 



We note that there are two quadrupolar contributions 
(with dimensionless moments Q'l, and Q^) which show 
an oscillatory behavior for radial distances r ^ h (i.e.. 



® In this respect the pictorial argument given in Ref. [13] 
(see Fig. 3(b) therein) is slightly misleading (at least in an 
asymptotic sense): there the assumed director configuration 
around the colloids for the case of a finite film thickness is a 
tilted dipole with a nonvanishing component v oc Pv^^x/z^. 
Actually, the broken symmetry in x-direction rather enters 
through a tilted quadrupole. 



zy 



^ V^™ 15 ^'^ {qorf 



Wi, 



(40) 



where J[Y]j+i/2{r) are the spherical Bessel functions of 
the first [second] kind and Yjm{d, (j)) are the usual spheri- 
cal harmonics for the standard set of spherical coordinates 
r, 9, (p. (The origin is again taken as the center of the cir- 
cular three-phase contact line.) The coefficients w ^^ are 
dimensionless multipole moments. 

The Dirichlet boundary conditions for v and w at the 
two interfaces can be fulfilled as before by constructing 
the full solution in terms of multipoles around the colloid 
and corresponding image multipoles as shown in Fig. [3l 
Since rotational covariance is broken by the parallel sub- 
strate anchoring, the solution for v contains a nonzero 
dipole contribution. Nevertheless, the director field still 
obeys a reflection symmetry with respect to the xz-plane: 
v{x, y, z) = v(x, — y, z) and w(x, y, z) = —w{x, —y, z). There- 
fore, the leading asymptotic behavior for v is given by 



For large film thicknesses the solution for the director field 
near the nematic-air interface in the regime i? ^ r <C /i 
should coincide with the solution for macroscopically thick 
films (Eq. (fTSJ) ). Near the nematic-air interface, one has 
V ~ Tix, w Ki Uy. Therefore one recovers the rotationally 



Q 



covariant quadrupole solution of Eq. ([TS]) for Q„ 
(obtained by equating Eq. p7|) with Eq. p5|) for i = 1, i.e., 
ri = x) and —?>^j2pnQ^/qQ = Q (obtained by equating 
Eq. dSZl) with Eq. ^ iov i = 2, i.e., ra = y). Since Q^ - 
qlQ and Q = 0{R^), Ql = 0{[qQRf) is a very small 
number. 

The Dirichlet boundary conditions for w at the sub- 
strate and at the nematic-air interface enforce that the 
contribution to w due to the quadrupole Q:^ and all cor- 
responding image quadrupoles is zero. This holds also for 
the contribution of all higher multipoles of degree j for 
which the radial dependence is given by Jj+i/2(^)- The 
only solution of the Helmholtz equation which fulfills the 
boundary conditions and which, as an additional require- 
ment, is smooth everywhere, is w = 0. Since the Bessel 
functions of the first kind are smooth everywhere, all re- 
spective multipole moments must be zero. (This does not 
hold for the multipole moments pertaining to the Bessel 
functions of the second kind since 1^+1/2 (?") is singular at 
r = 0.) 

The excess force on the system "colloid and nematic- 
air interface" follows from Eq. (j33p and can be expressed 
as 



2 



dA [v^^ — 2(7o w.z + w^^) 



(41) 



5sub 



utilizing the boundary conditions for the solutions v and 
w. For the multipoles appearing in the expansion of v 
the method described in Subsec. 12.2.11 may be used. In 
order to obtain the quadrupolar contributions to w we 
perform the summation over the image multipoles and the 
integration over the substrate surface numerically. Due to 
the slow decay of y5/2(?') ex (cos r)/r^/^ the stress integral 
is superficially divergent. However, a detailed asymptotic 
analysis yields convergence with the result 



AF ■ 



"^V^-ft7 



5C(5) Ql 
12 ft6 



+ c^{QlY\ ,(42) 



where Cw ~ 0.30 has been determined numerically. Note 
that the linear contribution due to the term ~ (jo v^z in 
the stress tensor (Eq. (|4T|) ) turns out to be zero, as can 
be easily checked by applying Gauss' theorem and the 
field equation Av = 0. Using the above considerations 
concerning the magnitude of the multipole moments P„, 
Qv and Q^ we obtain 



AF 



R 



4+2k 



R 



(43) 
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where ap and aq arc numerical coefficients of order 1 and 
K has been defined in Eq. ([38]) . Again the excess force 
on the colloid falls off rapidly with increasing film thick- 
ness such that the induced logarithmic capillary interac- 
tion (see Eq. ([H])) remains very weak. For the parameters 
characterizing the experimental system studied in Ref. [13| 
{R/h - 10"i, K/{"fR) - 10-'') we find Kicn(d) -- 
-[^g-ii-4K ^^y ln(_R/(i) which appears to be undetectably 
small. Note that for d < h the direct elastic repulsion re- 
mains essentially unchanged because in this regime the 
leading term of the elastic interaction is given by the re- 
pulsion between the two quadrupoles located at the colloid 
sites. In this case the image quadrupoles can be neglected. 



3 Discussion and conclusion 

We have investigated the effective potential between two 
colloidal microspheres of radius R floating at asymptoti- 
cally large distances d on an interface between a nematic 
film of thickness h and air (Fig.[T]). This effective potential 
is the sum of an elastic interaction caused by the director 
distortions around the colloids and of a capillary interac- 
tion mediated by surface deformations. We have analyzed 
the effective potential for large d employing the coarse- 
grained Frank free energy (within the one constant ap- 
proximation) for the director distortions and the linearized 
Young-Laplace equation for the interface distortions. 

In the case of a macroscopically thick nematic film, the 
director deformation around a single colloid is of quadrupo- 
lar type. Thus the induced elastic interaction between two 
colloids at distance d is repulsive and of quadrupolar type 
oc d^^. The capillary interaction is also repulsive and de- 
cays oc d~^ but it is much weaker than the elastic interac- 
tion. This rapid decay is a consequence of the mechanical 
isolation of the system, i.e., of the fact that the net force 
on the colloidal particles and the surrounding interface 
vanishes. 

A finite film thickness h leads to a violation of the me- 
chanical isolation of the system "colloid and interface". 
(The excess force is counteracted by the film substrate 
such that the whole experimental system is of course in 
mechanical equilibrium.) However, on the thermal energy 
scale the strength of the ensuing logarithmic capillary po- 
tential turns out to be very small. In the case of homeotropic 
boundary conditions on both sides of the nematic film the 
strength is proportional to (R/h)^"^ (see Eq. ([28|) and the 
subsequent discussion). For twisted boundary conditions 
(parallel at the bottom substrate, perpendicular at the 
upper nematic-air interface) there is a qualitatively dif- 
ferent asymptotic behavior of the director field due to a 
dipole contribution, which is induced by the broken rota- 
tional symmetry and vanishes for h —^ oo. However, even 
for these boundary conditions the strength of the loga- 
rithmically varying capillary potential remains extremely 
small, vanishing at least oc (R/h)^ (see Eq. (P5)) and the 
subsequent discussion) . 

Thus our analysis based on the mechanical isolation of 
the experimental system under consideration rules out a 
significant attractive contribution ^^ In d of capillary type 



to the effective potential between two colloids at a nematic 
interface. The amplitude of such a logarithmic contribu- 
tion vanishes rapidly for large film thickness h. Therefore, 
this effective pair potential does not provide a mechanism 
for the stability of isolated colloid clusters at a nematic 
interface as reported in Ref. [13]. 

Naturally one strives for other explanations for the ob- 
servation of stable clusters reported in Ref. [T3]. There 
the mutual center-to-center distance between neighbor- 
ing colloids in the cluster has been found to be between 
3R and 5R, depending on the radius R of the colloids. If 
the cluster stability is attributed to a minimum in the ef- 
fective pair potential at this range 3R . . . 5R of distances, 
the applicability of the above asymptotic considerations 
at such distances is doubtful for both the elastic and cap- 
illary contributions to the effective pair potential: 

— Elastic part: For the single-colloid problem, the asymp- 
totic, quadrupolar form of the director field (Eq. (fTS])) 
is based on the assumption that the deviations of the 
director from the preferred direction e^ are small: 
I'T^ajfyjl ^ 1. At a radial distance p from the center 
of the colloid this implies Q/p^ <C 1, which seems 
to be fulfilled for the dimensional estimate for the 
quadrupole moment Q = 0{R'^) and for the distances 
p — 3R . . . 5R under discussion. However, the absolute 
magnitude of the director deformations is not fixed by 
these arguments, so that it might be worthwhile to de- 
termine the minimum of the Frank free energy by a full 
numerical calculation in the presence of one or two col- 
loids at the nematic-air interface, similar to Ref. |38| 
where the explicit director solution around a single col- 
loid in the bulk has been determined. Our analysis in- 
dicates that only if the multipolar expansion of the 
director field around a single trapped colloid fails for 
distances p = 3R . . . 5R, there is a chance for attrac- 
tive elastic interactions between two colloids at these 
distances. Failure of the multipole expansion at these 
intermediate distances would point to a strong devia- 
tion of the magnitude of the multipole moments M„ 
of order n from the dimensional estimate M„ ~ i?"+^. 
Thus it is certainly worthwhile to determine multipole 
renormalizations for bulk and interface configurations, 
taking into account the full nonlinearity of the director 
field equations. 

Capillary part: Whenever the colloid-induced nematic 
stress deviates from the asymptotic, quadrupolar form 
(Eq. pT|l ). corrections to the asymptotic capillary po- 
tential will arise. However, due to the smallness of the 
dimensionless force ep (see Eq. (j9]), Ef ex K/{jR) « 
lO"'^) the perturbative treatment of the ensuing cap- 
illary deformations is justified. Therefore, one of our 
main findings, i.e., that the magnitude of the capillary 
potential (oc e|.) is always smaller than the magnitude 
of the elastic interaction (oc |£f|), is likely to hold also 
for elastic stresses deviating from the asymptotic limit. 
Another possible, subtle effect related to the anisotropy 
inherent to nematic phases allows the substrate to ex- 
ert also a torque on the colloidal particles and the sur- 
rounding interface. This case can be studied similarly 
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Pcap and the interaction en- 
,|/d)2. Thisde- 



to what we have done here. Actually, the result can 
be obtained straightforwardly if the analogy of cap- 
illary deformation with 2D electrostatics is employed 
(see, e.g., Ref. [2S|): The net force creates "capillary 
monopoles" Qcap with an interaction energy Vmon(c?) = 
Q2^p/(27r7) ln(CA/d) (see Eq. dH])). A net torque cre- 
ates "capillary dipoles" 
ergy between them is Vmcnid) ^ (|Pcapi 
cay is sufficiently slow to dominate asymptotically the 
nematic-mediated repulsion. However, dimensional anal- 
ysis yields Pcap ~ R Qcap so that the amplitude of this 
capillary energy, which must vanish for large film thick- 
nesses h —f oo, has also a too small numerical value to 
explain the experimental results reported in Ref. |13| . 

We recall that for simplicity our calculations have been 
carried out in the "strong anchoring" limit. "Weak anchor- 
ing" would lead to a different director field, and thus to 
different values of the dimensionless forces ep and e^^ (see 
Eq. ([5])). It would also lead to an explicit contribution 
to the expression (|10[) for the free energy accounting for 
the anchoring "wetting" energies (Eq. ^), which in the 
"strong anchoring" limit considered so far were subsumed 
in a (quantitatively negligible) additive renormalization of 
the surface tensions. The consequences of this more gen- 
eral case of no strong anchoring have to be explored yet, 
but we do not expect that this alters our conclusion that a 
logarithmically varying capillary attraction is ruled out by 
mechanical isolation (ep = Eu-), because this latter result 
is based on very general principles. 

Finally we note that for a two-dimensional system, the 
interactions between two upright circular colloids trapped 
at the nematic-isotropic interface have been studied nu- 
merically [39j using a Landau-de Gennes free energy. In 
three dimensions, this corresponds to long, cylindrical col- 
loids which are aligned side-to-side at the interface. For 
boundary conditions that yield a similar, quadrupolar be- 
havior of the director field around a single colloid [JD] , the 
effective interaction is found to be repulsive and consistent 
with a power-law decay (for intermediate distance d/R up 
to 7) . In this particular two-dimensional situation, the re- 
pulsive interactions appear to be longer-ranged (x d~^ and 
in the numerical results for the director field there is no 
trace of a sizeable interface deformation which would lead 
to capillary attractions. 

In summary, we have presented several arguments which 
rule out an asymptotic attraction of capillary origin in the 
effective interaction potential between two colloids at a ne- 
matic interface. If one discards effective pair potentials as 
the source of the stability of colloidal mesostructures as 
the ones found in Ref. |13| . the question arises whether 
this stability is a consequence of genuine many-body ef- 
fects [H]. Nematic surfaces (without colloids) can stabi- 
lize regular patterns of surface defects [30] . If these defects 
persist in the presence of colloids, the experimentally ob- 
served regular order of the colloids might be attributed to 
them. This issue calls for further experimental and theo- 
retical investigations. 
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A Corrections to the strong anchoring limit 

Here we consider the asymptotic director field around a 
single colloid at the nematic-air interface in the case of 
finite anchoring strength. For large distances from the 
colloid, the overall deviations from the preferred director 
n = (ni,n2,n3) = (0,0,1) are small and the total free 
energy of the nematic phase is given by the harmonic ap- 
proximation for the bulk part (Eq. (fT2|) ) and the Poulini 
expression for the surface part: 



^„c - f / d^r Y, iVn,f + -^ / dA{n- ba)' 

"'^no i = l,2 JAair-ne 

= f / d^r ((Vni)V (Vn2)') + 

^/ dAil-nj-nl). (44) 

The second equation holds because the normal e^ of the 
nematic-air interface ^air-no is parallel to the z-direction 
and n^ = 1. For particles with radii R in the ^m range, 
such as those investigated in Ref. [13], the dimensionless 
parameter a = K/{\Wi\R) is smaller than 1. Minimizing 
the free energy renders the Laplace equation in the bulk, 

An,^0 (i = l,2), (45) 

supplemented with the Robin boundary condition: 

d 

aR -77-ni{x, y,z = 0) ~ ni{x, y, z = 0) = (i = 1, 2) . 
az 

(46) 

The general, asymptotic solution is given by the multipole 
ansatz (Eq. P^ ) 



3 

E 

/3=1 



P, 






/3,7=1 



rp: 



(47) 



with the requirement of rotational covariance around the 
z-axis. In view of the boundary condition (Eq. (j46|) ) we 
note that the derivative in Eq. ([^S]) applied to a mul- 
tipole of order m generates a multipole term of order 
m -\- \. Thus, to leading order in 1/?', the boundary con- 
dition can only be met if the leading multipole (n^) ful- 
fills the strong anchoring condition n^(z = 0) = and 
is accompanied by a subleading multipole (rif^) which is 
connected to the leading multipole through the condition 
aRdzn\{z = 0) = nf{z = 0). Both the monopole and 
the dipole do not fulfill the strong anchoring condition 
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Fig. 4. A general configuration of the particle and the fiuid in- 
terface. The arrows indicate the orientation of the correspond- 
ing surfaces: Smcn is the interface, S A^J^ is the surface of the 
particle in contact with the lower (upper) fluid phase, S'waii is 
the surface of the container of the system (sketched here as a 
quadrangular box for simplicity) . The three-phase contact line 
between the particle and the interface is denoted as Co- 



n{{z = 0) = 0. Therefore, as before, we find that the mul- 
tipole expansion starts with the quadrupole solution and 
that it is necessarily accompanied by a hexapole: 






n. 



n 



s£ 






(48) 
(49) 



(50) 

Since 5^nf (z = 0) = 0, this solution (Eqs. (11H])-(ISD1)) ex- 
actly satisfies the boundary condition in Eq. ((46|) . Thus 
the asymptotically dominant director field consists of a 
quadrupole-hexapole superposition. Furthermore, the mag- 
nitude of the accompanying hexapole moment is small 
due to the factor a, therefore the results reported be- 
fore for strong anchoring (based on the leading quadrupole 
only) are unaffected by a finite anchoring strength at the 
nematic-air interface. 



B Force balance in a general configuration 

For the benefit of the reader, we discuss in this Appendix 
some previous results PSII^ concerning the force balance 
of a general equilibrium configuration and demonstrate 
how the amplitude of an asymptotic, logarithmically vary- 
ing interfacial deformation is determined solely by this 
mechanical condition of force balance. 

Figure [4] represents a colloidal particle in equilibrium 
at the deformed interface; in general the deformation is 
not small. The condition of mechanical equilibrium implies 



that locally the net force on any part of the system must 
vanish. Thus one has: 

1. Each of the fluid phases is in equilibrium. We introduce 
the forces exerted by each fiuid phase on the particle, 
on the whole fiuid meniscus, and on the wall of the 
container as 



ri(2) ._ 

part • 



ri(2) ._ 



^1(2) 
■ wall 



dA.ni(2), 






dA.ni(2), 



1(2) 



(51) 
(52) 
(53) 



respectively, in terms of the stress tensor n^^^)(r) in 
each fiuid phase with due account for the orientation of 
the surfaces (see Fig. 01). The superscript -'^^^^ indicates 
the lower (upper) phase in Fig.[4l The condition of me- 
chanical equilibrium of each phase under the infiuence 
of these three forces reads 



i.l(2) 
part 



ri(2) 



^1(2) 
■ wall 



0. 



(54) 



The total force exerted by the fiuids on the particle 



is F 



part 



part 



+ 



p2 

■^ part' 



and on the meniscus it is 
Fmon := F^(,„-|-F^jg„. The expressions for these forces 
reduce to those given in Eqs. ([6|) and ^, respectively, 
upon evaluation in the reference configuration depicted 
in Fig. H 

If the condition of mechanical equilibrium is applied 
locally to an infinitesimal volume in the bulk of each 
of the fiuid phases, it turns into V • n^^^^ = 0. In the 
nematic phase this condition yields the field equations 
determining the director field. 

The particle is in mechanical equilibrium under the 
combined action of Fpart and the tension exerted by 
the interface on the particle at the three-phase contact 
line Cq. This tension can be expressed in terms of a line 
integral involving the surface tension 7: 



If d£ec, 

'Co 



(55) 



where ec is the unit vector tangent to the interface, 
normal to the contact line, and oriented towards the 
particle side. Therefore, the condition of mechanical 
equilibrium reads 



part 



0. 



(56) 



Any piece 5'int C Saicn of the fiuid interface is in me- 
chanical equilibrium. The force on S'lnt exerted by the 
fiuid phases and the tension exerted on this piece at 
its boundary, dS-mt, are balanced: 






(57) 



with Gc oriented towards the exterior of S'lnt- If this ex- 
pression is applied locally to an infinitesimal piece of 
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interface, it turns into an equation for the interfacial 
deformation relating the mean curvature of the inter- 
face to the pressure jump accross it. If, in addition, 
the interface deviates only slightly from a flat inter- 
face (identified with the plane z = 0), this equation 
reduces in turn to the well-known equation ([59]) below 
for the local height u{x,y) over the plane z = 0. 

The net force balance of the whole system follows from 
the three separate balance conditions ([54|) . ([56|) . and ([57|) : 

with Sint = S'mon (sO that dSint = CqU Cwali; C'wall IS 

the three-phase contact line between phase 1, phase 2, 
and the container enclosing the system), one finds that, 
as expected, in equilibrium the net force of the outer en- 
vironment on the system must vanislu: 



^wall 



^ wall 



■7 



deer 



0, 



(58) 



where ec points to the exterior of the system. 

The condition of mechanical equilibrium can be ap- 
plied advantageously to obtain a precise statement about 
the amplitude of an interfacial deformation u[p,(j)) vary- 
ing logarithmically with the lateral distance p ^ R from 
the particle with radius R. Far away from the particle, in- 
terface deformations and their gradients are small, so that 
the linearized equation holds. Thus there is a distance ^ 
beyond which the linear theory is applicable: 



jVMr) = nUv) - nlir) (re^ext) 



(59) 



where the piece of interface 5ext is enclosed by the circle 
p = £, and the line Cwaii- The general solution to this 
inhomogeneous Laplace equation can be written as 



u(p,(/>) = Aa + Boln- + 
P 



where the last, approximate equality involves the leading 
order term in Vw of the line integral. (This approximation 
is justified because the interface deviates only slightly from 
a flat interface at the circle p = £.) Evaluating this latter 
term for the general solution (pDj) one finally finds that 
the amplitude of the logarithmic term in Eq. ((60l) is pro- 
portional to the force exerted by the upper and the lower 
fluid on the particle and on the meniscus: 



+ 00 

E 



-^m. T 



771—1 

1 

27^7 J So: 



B„ 



(60) 

COS7«(0 — 0„i) + 



h'dp' p'[nlJv')-nUv')] In- 



Bo 



1 
( 

27r7 



[^ part I ^ mcnj 



(62) 



In obtaining Eq. ([62]) we have used that in the region 5'oxt 
(where the interface is almost flat) one has to leading order 



dA-[U'^-U^]-e^ 



P 



d(t>dpp[nl{v)~nUr)]. 

(63) 
In conclusion, if the stress nl^{r) — Uzzi'"^) decays suf- 
ficiently fast with the distance p from the particle, the 
solution in Eq. (pD|) will be dominated asymptotically by 
the logarithm with an amplitude given by Eq. (|62p . pro- 
vided this amplitude does not vanish |3. As one can infer 
from Eq. (|54[) . this latter condition means physically that 
the walls exert a non-vanishing force: a logarithmic term 
can only arise if mechanical isolation of the system "col- 
loid -I- interface" is violated. (Note that Eq. ((62|) follows 
actually from Eq. ([58]) if the interfacial deformations are 
small.) 

We emphasize the generality of this result: it only re- 
quires that the interface departs slightly from a flat one for 
distances sufficiently far from the particle. (Otherwise it 
is actually not useful to speak of a logarithmically varying 
deformation to begin with.) In particular, the interfacial 
deformation close to the particle may be arbitrarily large, 
the particle itself need not be perfectly spherical, or it 
may even consist of a many-body conflguration lacking 
any kind of symmetry. 
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